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We clarify that the antiferromagnetic order in the distorted tetrahedral quasi-one- 
dimensional spin system induces electric polarization. In this system, the effects of low 
dimensionality and magnetic frustration are expected to appear simultaneously. We 
obtain the magnetic-field-temperature phase diagram of CU3M02O9 by studying the 
dielectric constant and spontaneous electric polarization. Around the tricritical point 
at 10 T and 8 K, the change in the direction of electric polarization causes a colossal 
magnetocapacitance. We calculate the charge redistribution in a small spin cluster con- 
sisting of two magnetic tetrahedra to demonstrate the electric polarization induced by 
the antiferromagnetism. 
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Since the discovery of a strong magnetoelectric effect in TbMnOs, 1 ^ multiferroics 
in transition-metal oxides have been extensively studied. 2 ) In the cases of the in- 
verse Dzyaloshinskii-Moriya interaction in a spiral spin structure 3-5 ) and the inverse 
Kanamori-Goodenough interaction in a collinear one, 6 ) the formation of the magnetic 
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superlattice plays an essential role in the magnetic-order-induced multiferroics. Geo- 
metrical magnetic frustration also plays an important role as the origin of a nontrivial 
spin configuration that breaks the spatial inversion symmetry. In this study, we demon- 
strate that the distorted tetrahedral spin system has the potential to show multiferroic 
behavior without any magnetic superlattice formation. We focus on the dielectric prop- 
erties induced by an antiferromagnetic (AFM) spin order in CU3M02O9 and discuss the 
possibility of the multiferroic behavior in an AFM spin system on the basis of the theory 
of the charge redistribution in frustrated Mott insulators. 7 ' 8 ) 

CU3M02O9 has two distorted tetrahedral quasi-one-dimensional quantum spin sys- 
tems of S — 1/2 spins along the 6-axis in its orthorhombic unit cell [Figs. 1(a) and 1(b)]. 
This compound has geometrical magnetic frustrations due to the tetrahedral spin align- 
ment and quasi-one-dimensionality simultaneously. This compound undergoes an AFM 
phase transition at Tn = 7.9 K without a magnetic field. 9 ' ^ Inelastic neutron scatter- 
ing measurements clarify the hybridization effects due to the J\ and J 2 superexchange 
interactions between two elemental magnetic excitations, i.e., that of the quasi-one- 
dimensional AFM spin system originating from the J 4 (= 4.0 meV) superexchange 




Fig. 1. Schematics of the distorted tetrahedral chain in CU3M02O9 along the &-axis (a) and in the 
ac- plane (b). The circles indicate the 5 = 1/2 Cu 2+ ions and the symbols distinguish their coordinates 
along the 6-axis. 2 ~ and Mo 4+ ions are omitted. The dashed, solid, bold, and dot-dashed lines 
distinguish the superexchange interactions J1-J4 between Cu 2+ ions. The solid rectangle in (b) denotes 
the unit cell, which contains two tetrahedral chains. 
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interactions and that of the isolated AFM spin dimers originating from the J 3 (= 5.8 
meV) ones. 11 ' 12 ) 

The AFM spin-wave branch arises linearly from zero energy at the magnetic zone 
center of (h, k,l) = (0, 1, 1), which is also the nuclear zone center below and above T^. 
This magnetic branch is strongly dispersive along the (0, k, 0) direction. This dispersion 
curve and the two-spinon continuum above mass gap energy A = 1.2 meV are typical 
features of a quasi-one-dimensional AFM spin system. There are no indications of the 
formation of commensurate or incommensurate magnetic superstructures, at least at 
our experimental resolution. 

To obtain the magnetic-field-temperature (H-T) phase diagram, we measured the 
T and H dependences of the dielectric constant e a and the electric polarization P a 
under an electric field E a along the a-axis (a = a or c) of CU3M02O9. In this letter, 
we focus on the effects of the magnetic field along the c-axis. We prepared a platelike 
single crystal of CU3M02O9 whose cross section and thickness are typically about 60 
mm 2 and 0.4 mm, respectively. To form a capacitor, the faces were coated with gold and 
attached using two gold wires. The capacitance, of which the typical value was on the 
order of 10 pF, was measured using an impedance analyzer (Yokogawa-Hewlett-Packard 
4192A). e a was obtained from the capacitance at 100 kHz with a peak voltage of 1 V. 
To confirm the ferroelectric behavior, the electric polarization-electric field loop (P a -E a 
loop) at 1 Hz was recorded using a modified Sawyer- Tower circuit with a peak voltage of 
about 200 V. The magnetic field was applied using a superconducting magnet (Oxford 
Instruments, Teslatron S14/16), of which the maximum magnetic field was 16 T, and 
a variable-temperature insert cryostat was used to set the temperature from 1.5 to 300 
K. 

Figures 2(a) and 2(b) show the typical T dependences of e a {ot = a or c) under 
a fixed H (e a -T curves), respectively, each of which has a (local) maximum eP eak at 
jipeak r-j-^g va i ues f j-peak a g a i ns t H are plotted in the H-T phase diagram by the solid 
symbols in Fig. 3. Figures 2(c) and 2(d) show the typical H dependences of e a at 8 and 
4 K (e a -H curves), respectively. At 8 K, as shown in Fig. 2(c), the e a -H and e c -H curves 
have one and two peaks, respectively. These are plotted in the phase diagram by open 
symbols in Fig. 3. We observed a colossal magnetocapacitance (800%) at 8 K in Fig. 
2(c). Here the maximum e a slightly above 10 T was compared with the minimum e a at 
approximately 6 T. This magnetocapacitance is larger than the maximum one reported 
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Fig. 2. Typical temperature dependences of the dielectric constants under fixed magnetic fields [(a) 
and (b)] . For visibility, the data were shifted. The magnetic held dependences of the dielectric constants 
at 8 and 4 K are shown in (c) and (d), respectively. 



before. 13 ^ This is probably due to strong fluctuations around the tricritical point at (H, 
T) = (10 T, 8 K) in the H-T phase diagram. At 4 K, as shown in Fig. 2(d), e c gradually 
decreases with increasing H above 6 T. At about 8 T, a small peak is observed in the 
e c -H curve, which is indicated by an arrow and is plotted in Fig. 3 by an open symbol. 
Around this H, e a shows a sudden increase. 

At 4 K, as shown in the inset of Fig. 2(d), the e c -H curve under a zero- magnetic- 
field cooling process rapidly decreases with increasing H from T. This anomaly is not 
observed during the reduction of the magnetic field. We plot H, where this hysteresis 
effect disappears, in Fig. 3. Around this H, the H dependence of the magnetization M c 
(the M c -H curve) along the c-axis shows a small jump with a magnetic-field hysteresis 
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Fig. 3. i7-T phase diagram of C113M02O9. The shape of the symbols distinguishes the physical 
quantities to be used to obtain the phase boundary. The triangles, squares, and circles denote the 
dielectric constants along the a- and c-axes and specific heat, respectively. The solid (open) symbols 
denote the phase boundary obtained from the data of the T (H) dependence. 



effect. 9 ^ In some multiferroic materials, a change in electric polarization accompanied by 
a jump of the magnetization has been reported. 1 ' 14 ) These facts and the present result 
showing that the peaks in the e a -T curves become sharp at 6 T [Figs. 2(a) and 2(b)] 
suggest that the zero-field ground state contains some fluctuations, which is consistent 
with our picture describing the canted antiferromagnetism at zero and finite magnetic 
fields. 9 ' 15 ) In this picture, the canting direction of the spin moment contains randomness 
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Fig. 4. Typical polarization-electric-field loops. The two loops at 4 and 10 K under T in (a) almost 
overlap with each other. 



after a zero- field cooling process. The jump of the magnetization in finite magnetic fields 
is understood as the alignment of the canting direction. Once it occurs, it survives even 
at a zero magnetic field because of the internal magnetic field, indicating the possibility 
of the magnetic-field hysteresis effect in this system. 

Together with the additional data and phase boundary obtained from the T de- 
pendence of the specific heat under H, 10 ^ we obtain the H-T phase diagram in Fig. 3. 
One can see that the increase in T^ eak at 6 T [Figs. 2(a) and 2(b)] corresponds to the 
change in T N under H. 10 ^ We emphasize here that the peak in the e a -T (e a -H) curve is 
not necessarily in agreement with the critical temperature (magnetic field) of the phase 
transition. Therefore, it is natural that the plots in the H-T phase diagram obtained 
from the data of e a do not perfectly trace the phase boundary obtained using specific 
heat. We conclude that the phase diagram contains four different phases. 

To observe the spontaneous electric polarization directly, we measured the P a -E a 
loop. Figures 4(a) and 4(b) show the typical results at various T and H. At T, as 
shown in Fig. 4(b), the P c -E c loop below T N shows clear ferroelectric behavior. The 
spontaneous electric polarization density (~500 /iC/m 2 ) is close to that of TbMnOs. 1 ^ 
As shown in Fig. 4(a), the P a -E a loop at T does not show ferroelectric behavior. 
Thus, we conclude that CU3M02O9 at T below Tn is in the ferroelectric phase with 
spontaneous electric polarization along the c-axis [the FE(P S / /c) phase]. This result is 
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consistent with the following facts: ej? eak ~ 10 x e^ ca,k in Figs. 2(a) and 2(b) and a strong 
microwave absorption was observed at approximately T N at a zero magnetic field. 

Figure 4(b) shows the H dependence of the P c -E c loop. At 5 T, the P c -E c loop 
indicates the FE(P S / /c) phase because it is similar to that at T (data not shown). 
Instead of the closed P c -E c loop at 13 T, the ferroelectric P a -E a loop was observed at 
14 T. The eP eak at 13 T becomes about ten times larger than fP cak at 6 T [Fig. 2(a)]. 
The eP cak at 13 T becomes about ten times smaller than eP eak at 6 T [Fig. 2(b)]. Judging 
from these results, we confirm the FE(P s //a) phase at 13 T. The spontaneous electric 
polarization density (~800 /xC/m 2 ) is comparable to that of DyMn0 3 at 21 K. 13 ) 

As shown in Fig. 2(d), both the e a -H and e c -H curves show anomalies at approxi- 
mately 8 T at 4 K. As shown in Fig. 4(a), the P a -E a curve suggests a double hysteresis 
loop at 7.9 T. This effect is explained by the electric-field-induced electric polarization 
and indicates strong fluctuations around the phase boundary. We conclude that the 
change in the direction of the spontaneous electric polarization occurs at the phase 
boundary running from (H, T) = (8 T, 2 K) to (10 T, 8 K). At approximately 8 T, 
a change in the electron-spin- resonance spectrum was reported. Together with the 
phase transition temperatures obtained from the T dependences of M c under various 
if, 10 ) we conclude that this compound is a multiferroic material. 

In many cases, the origin of multiferroicity has been discussed on the basis of the 
formation of the magnetic superlattice. Unfortunately, the spin structure below Tn 
at T has not been clarified yet. However, we emphasize again that the AFM long- 
range order below T N has been established by the magnetic dispersion curve obtained 
from inelastic neutron scattering. 12 ) Moreover, the anisotropic magnetization has been 
quantitatively explained on the basis of a weakly canted AFM long-range order. 9 ) 

The point of discussion is how the AFM long-range order induces ferroelectricity. 
In the following, we focus on the charge redistribution effects caused by the three- 
spin ring exchange interaction in geometrically frustrated Mott insulators, proposed by 
Bulaevskii et al. 7 ' 8 ^ We demonstrate that the spin cluster in Fig. 5 has the potential 
for multiferroic behavior as a result of the antiferroelectric (AFE) polarization at zero 
electric field and the ferrielectric (FRE) polarization under an electric field. This is the 
minimum system that maintains the tetrahedral spin arrangement and the inversion 
center at spin site #4 simultaneously. The formation of the magnetic superlattice is 
impossible. 
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Fig. 5. Six-spin cluster under periodic-boundary condition. indicates the ith 5=1/2 spin. This 
spin cluster can be divided into two tetrahedral spin clusters, each of which has an electric polarization 
Pi or P 2 . The charge redistribution at the center site (#4) is divided so that both Pi and Pi exist 
on the bisector between sites #1 and #4. 

We treat this system as a Mott insulator. The electronic band from d x 2_ y 2 orbitals 
in Cu 2+ ions is half-filled. In the limit of the strong on-site Coulomb repulsion U, where 
the charge degrees of freedom are frozen, the system can be mapped on the quantum 
spin system. The spin Hamiltonian is given by 



where the sum in the first term on the right-hand side runs over all the possible spin 
pairs (Si and Sj at the #i and #j sites, respectively) connected through the hopping 
parameter tij. The values of 4i?-/i7 correspond to the exchange interactions in Fig. 5. 
The second term on the right-hand side is the magnetic energy from the local mag- 
netic field H x ° c at the #i site, where g and //b are the g factor and Bohr magneton, 
respectively. 

On the basis of the interchain interaction in CU3M02O9, 12 ) we set the values J\ = 
J 2 = 1 meV, J 3 = J3 = 5.8 meV, J 4 = 4 meV, and g = 2. We introduced the staggered 
magnetic field Hafm = |-H"!° C | along the quantization axis working on only the #1 
and $4 sites. Hafm induces the AFM spin order at the site on the spatial inversion 
symmetry, i.e., this site is on the anti-inversion center in the term of the magnetic space 
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Fig. 6. Spin moments projected on quantum axis (a) and charge-redistribution parameters (b) in six- 
spin cluster under periodic-boundary condition as functions of staggered magnetic field. The parameter 
sets that give the antifcrroelectric (AFE) and ferrielectric (FRE) phases are given in the text. The 
electric polarizations are shown by the arrows. 



group. Using the exact diagonalization, we calculated (Si) and the charge redistribution 
5rii at spin site #?, 7 ' 8 - ) 

5n t = ]T + S k ) - 2S r S k ] , (2) 

(i,j,k) 

in the periodic-boundary six-spin cluster at K. Here, the sum on the right-hand side 
runs over all the possible spin triangles connected through exchange interactions. The 
amplitude of the electric polarization |Pi(2)| is proportional to Sni^) — 5n 2 (i) — 5n^) + 
5n 4 (6), 

As shown in Figs. 6(a) and 6(b), |(<S'{i,4})z| and |P{i,2}| increase larger with in- 
creasing Hxfm- Even in this case, the spins at the sites #2, #3, #5, and #6 still form 
nonmagnetic spin dimers. We plot 5nU" 3 ^ 2 in Fig. 6(b) because we could not precisely 
estimate U, and calculations in larger systems are necessary to discuss the origin of the 
finite 1^(1,2}! without the staggered magnetic field because of the strong system-size 
dependence. 

Under a uniform electric field, the exchange interaction slightly changes as a result 
of the breaking of the spatial inversion symmetry. 17 ) If the exchange interactions of the 
nonmagnetic spin dimers are different (J3 = 5.0 meV and J3 = 6.6 meV), we obtain a 
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finite net electric dipole moment as a result of the FRE alignment of the electric dipoles 
(Pi 7^ — P 2 ) even though (Si) is not changed. From this calculation, we conclude that 
the tetrahedral quasi-one-dimensional spin system has the potential to exhibit AFM- 
AFE and AFM-FRE type multiferroic behaviors. The finite net electric dipole moment 
appears when the time and spatial inversion symmetries are broken. 

In this study, we demonstrate the possibility of electric polarization induced by the 
AFM long-range order in a distorted tetrahedral quasi-one-dimensional spin system. As 
experimental evidence, the H-T phase diagram of CU3M02O9 under a magnetic field 
along the c-axis was shown by studying the temperature and magnetic-field dependences 
of the dielectric constant. The ferroelectric behavior was observed in the polarization- 
electric-field loop and a change in the polarization direction was observed. Around 
the tricritical point at (10 T, 8 K), a colossal magnetocapacitance was observed. We 
showed that a six-spin cluster that corresponds to a tetragonal quasi-one-dimensional 
spin system has the potential to become an AFM-AFE multiferroic state without the 
formation of a magnetic superlattice. The possibility of the AFM-FRE multiferroic state 
under an electric field was discussed. 
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